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Abstract 


In order to represent the magnetic field B in the magn et 0 - 
sphere various auxiliary functions can be used: the current den- 
sity j. , the scalar potential j , toroidal and poloidal poten- 
tials '4 / 1 and and Euler potentials (oi, j&) — or else, 

the components of B may be expanded directly, with constraints 
ensuring the vanishing of V»B . The most versatile among the 
linear representations is the one based on ( y , S' ) ; it has 
seen relatively little use in the past but appears to be the most 
promising one for future work, Euler potentials are non-linear 
and can only be recommended for cases where their special pro- 
perties are utilized, e.g. the representation of electric poten- 
tials when En = 0 o Other classifications of models include 
simple "testbed" models vs. "comprehensive" ones and analytical 
vs. numerical representations. The electric field E in the 
magnetosphere is generally assumed to vary only slowly and to be 
orthogonal to B, allowing the use of a scalar potential 
which may be deduced from observations in the ionosphere, from 
the shape of the plasmapause or (as Mcllwain has done) from 
particle observations in synchronous orbit. A simple model of 
<f> is discussed and general implications are described. 



This talk is meant to he a review of technical points - of methods 
and ideas - involved in the construction of quantitative models of 
magnetic and electric fields in the magnetosphere. 

Because time is limited, I shall not devote ny talk to the cataloging 
and comparison of existing methods: I have a review article available 
which does just that for models of the magnetic field and you are welcome 
to take a copy with you, to read on the plane home. There are also 
available some copies of a somewhat more restricted piece of work on 
electric fields. 

Instead, I would like to use the time to bring a bit of order to 
the profusion of models - to classify the wide variety of models accor- 
ding to mathematical type , representation and application . When we use 
a model our choice generally depends on the application for which it is 
intended and this classification, I hope, will make it clearer what is 
available. 

(Figure l) 

Let me start with the magnetic field . One basic 
classification depends on the auxiliary functions which are used for 
representing the field. 

The first two of the representations shown are based on 
the current density j , which is generally introduced in one of two 
ways. First , there exist cases where is an observed quantity - 
say, the tail's current sheet, field aligned currents or the ring current 
as deduced from particle populations. Akasofu and Chapman, for instance, 
carried out extensive work on ring current fields based on this approach. 
Of course, what you get is then a model of what your theory predicts the 
field to be, not necessarily a representation of B as observed. 

Secondly , you can express in some general way and fit 
the expansion coefficients so that the observed field B is represented 
as closely as possible: this is Bill Olson's approach. By using j one 
can assure that the divergence of B in the derived model vanishes, 
although at first sight this does not appear much of an advantage, since 
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j[ also has to be divergence-free . The real advantage is that while j 
may be a rather discontinuous current density distribution - current 
filaments or what not - the resulting B is rather smooth, since it 
is obtained by integration. 

Next, you have the representation by a scalar potentia 

7 : this is only good for curl-free fields, but it is the preferred 
method for representing the main geomagnetic field which originates in 
the earth's core. Usually 7 is expanded in spherical harmonics and the 
number of terms can go up to a 100 or so, depending on the accuracy 
which you want. 

The next method uses two functions V . and ^ which, as far as 

■ ■ - • 1 1 2 . 

I know, have no names, and I’ll therefore call them here t oroidal 
and poloidal potentials, respectively, since the 
terms in 'which they appear are called the toroidal and poloidal compo- 
nents of B * This very powerful representation - it is equivalent to the 
use of spherical vector harmonics - was introduced into dynamo theory by 
Walter Slsasser about 30 years ago and is well-known to astronomers, but 
not, apparently, to those engaged in magnet ospheric physics, It deserves 
more attention from us and 1 will have something to say about this later 

Olio 


To give you some intuitive feeling for what these functions mean, 
notice from the second line in item (4) that the toroidal field is 
perpendicular to r : it thus represents field lines circling the 
origin in some manner, like field lines which circle a wire in which an 
electric current flows , 

The poloidal component, on the other hand, resembles what you find in 
magnet ospheric models * The dipole field, for instance, is poloidal; more 
generally, you will note that if is harmonic, the poloidal component 

is curl-free and, in fact, all curl- free fields can be thus represented. 
The representation is unique - that is, any part of B has to be either 
poloidal or toroidal, there remains no ambiguity. 
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The preceding 4 representations all form a single group: they are 
ail linear and can therefore "be superposed as -we see fit: we could, for 
instance , combine a main field represented by y with a tail field 
given bjr j - this is done in the Mead-Williams model - and improve 
the fit by adding expansions of ^ and ^ ° By contrast the next 

method on the list - Euler potentials -is not linear, 
since in using it you multiply derivatives of 0^ by those of ^ . 

Because of this non-linearity one cannot in this case add up contri- 
butions - instead, c < and ^ have to be calculated from the beginning 
for the total field. This is a great inconvenience, so unless you have a 
very good reason - or work with the dipole field, where oL and are 
simple - it may be better to use a different representation. 

The advantage of Euler potentials over other methods 
is that they give an explicit analytical representation of magnetic 
field lines. Whenever the physics of the situation demands such a 
representation, they tend to be extremely useful: later on, when 
electric fields are discussed, we shall see one example of this. 

Finally, B can be expanded in a general analytical or numerical 

way. 

(Figure 2) 

One problem here is in ensuring the vanishing of V'B . The 
magnetospherie models of Mead and Fairfield, for instance, expand the 
components of B in powers of x, y and a, as shown in the figure, 
and they ensure the vanishing of V °B by the addition of linear 
constraints, which are taken into account (when the coefficients 
are derived) by the method of Lagrangian multipliers. 

Notice, however, that the same result could be obtained more neatly 
if we used toroidal and poloidal potentials and expanded them in 
powers of (x, y, z). If you do this, then ^ r and are sums 
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Direct Representations of B 

G«>H» Mead and D.Ho Fairfield, "A Quantitative Magnetospheric 
Model Derived from Spacecraft Magnetometer Data 1 ', JGR 80 , 523 
February 1975 > use the representation 


_ / km n-K-m 

B. = L , a.. x y z 

1 tTTTA lKmn 

jS ^ HI j J1 


n — 2 


Linear constraints assure the vanishing of V°B and 
terms are omitted to preserve symmetry. 


However, if 


B = V X r +- VXV,XH-'r 

r 2*— 


k ui n-k-m 



k,m,n 


a., x y z 
ikmn 


then, since 


A A A 

x X r y y f- z z 


the same expansion results vith no need for constraints, 


For "better control at large r it helps to modify 


the expansion to 


r - * 

V. =Z . a k~ y ffi z n-k-ra e -(r/r 0 ) 

1 k,m,n ikinn 


Figure 2 
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of the unit vectors in the (x, y, z) directions multiplied ty polynomials 
in (x, y, z), and if you take the curl or double curl you are still left 
with expressions of the same sort* 

This approach not only eliminates the need for constraint equations 
hut also makes it easy to generalize the method,, The Me ad- Fairfield 
expansion stops at quadratic terms - those with n = 2 - because higher 
powers are hard to control near the boundary and besides, the constraints 
become non-linear. However, with ^ x and 4> a you can add an expo- 
nential term which limits the expansion terms at large distances, giving 
a model similar to the one devised by Olson but with strict control over 

V °B . 

Some time ago I have developed a computer program which implements 
this method and it seems to work quite well. If anyone here is interes- 
ted, I will be glad to discuss it in private later on. Other modifications 
to this approach could also be devised: because toroidal and poloidal 
potentials are such versatile tools, I expect them to be important in 
future development of quantitative magnet ospheric models. 

(Figure 5 ) 

With so many methods of representation available, many different 
models can be - and have been - constructed. They seem to fall into two 
main classes. There are “testbed models" which aim 
at simplicity: you use them in theoretical work when you want to investi- 
gate effects involving some qualitative properties of the field without 
dragging in too much complexity. For instance, if you wish to develop 
a theory of effects due to the South Atlantic anomaly, you might be 
satisfied - at least at first - with the eccentric dipole. 

On the other hand - and of more interest to this meeting - there 
exist "comprehensive models" (some people here 
may call them "quantitative models", although strictly speaking all 
models discussed here are quantitative) - which try to represent 
observations as accurately els possible. The procedure by which such 
models are derived usually involves some mathematical representation 



- 7 - 


fl>. 


Types of Geomagnetic Models 


Class 

of Models 

Specific 

Specific 

Models 

Applications 

i ... .... — - i 



Dipole field 

General - simplest 

T e s 

tied 


approximation 

M o 

dels 

Eccentric dipole 

South Atlantic anomaly 



Image dipole 

Simple model of 



Mead's 3-term 

distorted field 



model 

and dayside boundary 



Mead-Williams 




3 -term ( ^ ) 




2 - dimensional 
models of tail 
field 

Particle motion 
in plasma sheet 


Comprehensive 

Models 


Main field 7 


Study of internal 
field and of field near 
surface of earth 


Main field (o£,p) Conjugate points 


Olson and Pfitzer 
model 

Mead and Fairfield 
model 

Generalization by 
V 1 and 

Magnetospheric 

(<*, p> 

(in future?) 


General use* 

Fitting of satellite 
data. 

Correlation ^with 
tilt, Kp, sector etc. 


Mapping of E 


Figure 3 
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which contains a number of unknown coefficients , and the values of the 
coefficients which best fit the observed data are derived by least squares 
fitting. Such models perform several useful functions: 

(1) They average out fluctuations in the data, 

(2) They help relate observations of particles etc. to the "real" 

magnetosphere . 

( 3 ) They enable one to extract from large data sets the average behavior 

of the magnetosphere - how it changes with K p , with the tilt angle 

of the dipole axis, with the interplanetary field and its sectors, 
with solar wind pressure, and so forth. 

A word of caution is however appropriate: such models do not provide 
data where none is available. It is the nature of models to bridge over 
regions of sparse data, or to extend to distances beyond those for which 
data exists, and the model is then no more than a mathematical inter- 
polation or extrapolation. This is especially important to remember with 
models of the electric field, like Mcllwain's - and even if the author 
there warns all users that the model is only valid in a limited region, 
there exists great temptation to follow it beyond its limits. 

Time does not allow me to go into other details, but there exists one 
more division of models which should be discussed, namely of analytical 
vs. numerical . 

■In all representations given in figure 1 the functions representing 
the field may be given either analytically or numerically. So far almost 
all models have been analytical, simply because even a sparse numerical 
grid introduces a tremendous jump in the number of coefficients handled. 
Yet we might be approaching the limit of practical accuracy in analytical 
representations . 

One simple remedy is to use different representations for different 
regions and splice them together where they meet. This might be a useful 

The facing page is left blank so that 
figures will remain matched to the text 
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thing to do in the tail region - one really should not expect the same 
expansion to describe the tail and the main magnet osphe re , since the 
physical processes in the two regions are quite different. Ultimately, 

I suspect, there exists a stage at which it pays to express the difference 
between the average field and our "best" analytical model by a coarse 
numerical grid - since the difference would be small, the grid would not 
require great accuracy. At the present time, however, the dispersion of 
our observational data is so large that analytical models satisfy all 
our needs. 

Let me now switch over to the electric field E 
where things are in a much worse shape, mainly because of the lack of 
data. 

(Figure 4 ) 

I do not have the time here to go into the history of this 
subject, which is fascinating, or the theory, which is controversial 

- let me just say that AlfVen's original speculation about a large- 
scale dawn-to-dusk electric field across the magnetosphere seems to be 
borne out. The same electric field probably also extends 
across the geomagnetic tail, while near the earth it must be modified by 
the addition of an extra component due to the earth's rotation (and 
perhaps some contribution from ionospheric motions). 

In most applications it can be assumed that the magnetic field does 
not vary with time, so that E can be represented by a scalar poten- 
tial <f> , If conductivity along magnetic field lines is high then such 

lines will be electric equipotentials: the result is best expressed 
when the magnetic field is given in terms of Euler potentials (in fact, 

I know of no other way) and reduces to ^ being a function of and 

alone . 

For instance, the contribution of the co-rotation field to , in 
the case of an axisymmetrical model of the main magnetic field, is 



- 10 


Models of the Electric Field 


For time- independent fields 

E = - 


If E«B = 0 then 

^ = <t> (ot,^ ) 

and <f> in the magnetosphere is determined by its value in 
the ionosphere or (for closed field lines) in the equatorial 
plane o 

If B is symmetric around the axis of rotation, the 
contribution of co-rotation to $ is 

i - 0< CO R 

; T cor 0 e 

i where (*) is the angular velocity of the earth *s rotation 

! and R is the earth’s radius » 

| e 

I 


Figure 4 
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o4 tL>R e 


as shown in Figure 4o If the field T s asymmetry _is taken into account , - 
the rotation of the earth leads to a finite ^B/^&t in the frame of 
reference of the magnetosphere and one cannot use ^ alone any more® 
Ways do exist for handling this situation hut I do not have the time 
to describe them. 

If <p is expressed in terms of the field-line parameters q£ and 
one only has to know its value at one point on each field 
line in order that be fully specified# Convenient choices for that 
point are either at the "roots rr of the field line in the ionosphere 
or in the equatorial plane; as it turns out, these are also the two 
locations where most of the information about E is obtained* 


(Figure 5 ) 

The electric field in the upper ionosphere has been inferred from 
ionospheric currents, barium cloud drifts, auroral motions and direct 
observations from 0G0 6 and Injun 5* from rockets and even from balloons, 
and all the evidence points to a two-celled electric field as shown 
in Figure 5 * What the figure shows is a schematic map of equipotentials 
in the polar cap, and below it you can see a sketch of how the dawn-dusk 
component of E varies during a pass over the middle of the polar cap. 

If, in the map drawn here, one introduces plane polar coordinates 
(R, i* ), then ^ can be represented (very nearly) by the analytical 
functions given on the slide. These functions contain one adjustable 
parameter k which represents the steepness with which the electric 
field falls off just outside the polar cap boundary: from profiles of 
the polar electric field, similar to the one drawn in Figure 5 and 
obtained by Heppner on 0G0 6 , one finds that k « 4 * Note that it 
is the region outside the polar cap that interests us most, 
since it corresponds to field lines which close inside about 10 R e * 
Field lines connected to the polar cap are either open or lead into the 
tail and are much harder to include in a model, since their properties 
are not well known* 
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Figure 5 
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(Figure 6) 


You no™ translate your map into (ot*^)* add the co-rotation field 
and thus get a model valid for the entire volume threaded by your 
field lines * You can,, for instance* map the electric field into the 
equatorial plane and it is interesting to note what happens: with 
k = 2 the field there - without co-rotation - is a constant field 
from dawn to dusk* with equipotentials stretched along the noon- 
midnight direct ion o With k less than 2 the equipotentials are pinched 
near earth while with k more than 2 - the actual case - they bulge 
out there » The sketches at the bottom of the figure show how it all 
looks when co-rotation is added* 

The method outlined here is probably the most feasible for mapping 
out E in any detail: in 5 years or so* if the Elect rodynamic Explorer 
satellite ever becomes reality* we ought to be getting quite detailed 
maps of the electric field in the polar ionosphere as functions of 
( o(* ^)* and they can then be mapped into the equatorial plane or anywhere 
else 9 

The next figure ( Figure jO shows how the k = k 
equipotentials actually look in the equatorial plane* 

The closed contour marks the boundary at which the co- rotation field 
becomes dominant and this seems to correspond to the plasmapause* 

Volland (JGR ^8* 171* 1975) used the observed shape of the plasmapause 
and also obtained k ** 4 for quiet times* which seems to support 
this approach* For disturbed times he got k » 2*75 and in addition 
there was a slight rotation of the pattern* so that it was no longer 
symmetric with respect to the noon-midnight meridian* This rotation 
shifts the bulge of the plasmapause towards midnight* as has been ob- 
served* Earlier calculations of this kind* by Vasyliunas and by Nagata 
and Kokubun* are also cited by Volland* 
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In the equatorial plate E itself is too "weak to be measured 
directly* but Mellwain has constructed some rather detailed models of 
E based on observations by ATS-5 in synchronous orbit. His data 
come from enhanced fluxes of low energy particles - especially 
protons - injected during substorms* and he assumed that all particles 
were impulsively injected at the inner edge of the plasma sheet* at 
a single instant. He also assumed that the electric field did not vary 
in time and proceeded to express its potential by means of a general 
mathematical expansion: the coefficients of this expansion were adjusted 
until they fit as closely as was possible the observed particle spectra 
and the times at which they were observed. 

This method claims to give E within the range of 5 to 10 Rg * 
although it is difficult to assess its accuracy. I hope that later in 
the session we will have the opportunity to hear more about it. 

Ultimately* in models of both B and E * we are going to run 
against the limit imposed by the variability of these fields. The vari- 
ation of E is especially pronounced and has been explored by Chen* 
Grebowsky and others: they deduce it from variations of the plasma- 
pause, which lead, to "tails" and/or "islands" of plasma isolated from 
the main body of the plasmas phe re . 

There still remains a lot here that’s not only poorly mapped but 
also poorly understood. I hope that within the next 5 years we will 
obtain at least good models of the average magnetic and electric fields 
in the magnetosphere: after that we might try what happens on shorter 
scales of space and time. 


Thank you 



